Local Measurements of Magnetization in Mni 2 Crystals 
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The spatial profile of the magnetization of Mni2 crystals in a swept magnetic field applied along the 
easy axis is determined from measurements of the local magnetic induction along the sample surface 
using an array of Hall sensors. We find that the magnetization is not uniform inside the sample, 
but rather shows some spatial oscillations which become more prominent around the resonance field 
values. Moreover, it appears that different regions of the sample are at resonance at different values 
of the applied field and that the sweep rate of the internal magnetic induction is spatially non- 
uniform. We present a model which describes the evolution of the non-uniformities as a function 
of the applied field. Finally we show that the degree of non-uniformity can be manipulated by 
sweeping the magnetic field back and forth through part of the resonance. 



Molecular magnets, or single molecule magnets, 
arc typically composed of magnetic cores sur- 
rounded by organic complexes. With a total spin 
of S = 10 and strong easy axis anisotropjii* 2 * 3 ^*^^, 
[Mn 12 Oi2(CH 3 COO)i 6 (H 2 0)4] • 2CH 3 COOH • 4H 2 
(generally referred to as Mni2-acetate) has received 
a great deal of attention. When crystallized, Mni 2 
forms a tetragonal lattice^ with the easy magnetization 
direction along the c-axis. The magnetic cores are well 
separated^, so that exchange interactions between the 
molecules are negligible^LiS. The main terms in the 
spin Hamiltonian of Mni 2 are given by 



H = -DSi 
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The first term is the magnetic anisotropy energy, with 
D = 0.548(3)K while the third term represents the 
next higher-order term in longitudinal anisotropy, with 
A = 1.173(4) x lQ-SR 11 ! 12 ! 13 . The second term is the 
Zeeman coupling to a magnetic field B z along the easy- 
axis. In the absence of a magnetic field, this Hamilto- 
nian can be represented as a symmetric double well po- 
tential with a set of energy levels corresponding to the 
(2S + 1) = 21 allowed values of the quantum number S z . 
A magnetic field B z lifts the degeneracy of ±5 Z states 
on opposite sides of the potential barrier, and a strong 
field B z causes most of the molecules to occupy states in 
one of the wells. When B z is swept, e.g., from positive 
to negative values, molecules relax from one side of the 
potential barrier to the other either via thermal activa- 
tion over the barrier or by quantum tunneling through 
the barrier. Measurements below the blocking tempera- 
ture of 3K show resonant tunneling of the magnetization, 
manifested as a series of steep steps in the hysteresis loop 
at roughly equal intervals of magnetic fieldi2ii£*i£ii£. The 
steps occur at magnetic field values where states on op- 
posite sides of the barrier have the same energy, and are 
at "resonance". 

One of the common experimental techniques to study 
the tunneling of the magnetization is to measure the mag- 
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FIG. 1: Schematic diagram of the experimental setup. The 
Mni2 crystal, with a typical size of z = 400 )im, y — 100 \im 
and x = 40 /im is placed directly on the surface of an array 
of 10 x 10 fim 2 sensors that are 10 /jm apart. The sensors 
cover only half the sample and measure the perpendicular 
component (B x ) of the magnetic induction at the surface. 



netic response of the crystal when an external magnetic 
field H a is applied parallel to the easy axis, and is swept 
through a series of resonances along the hysteresis loop. 
To date, studies of extended samples have provided the 
magnetic induction (and its variation in a swept magnetic 
field) as averages over the whole sample, or they have as- 
sumed that values measured in a limited region of space 
represent the average. As such, these measurements do 
not provide any information on the spatial profile of the 
magnetization in the sample, or on how the relaxation 
process propagates spatially within the sample. 

The experiment reported here utilizes an array of sen- 
sors to measure the spatial variation of the magneti- 
zation in a Mni2 crystal as H a is swept. Within our 
model we find the magnetization to be significantly non- 
uniform along the sample, with the non-uniformity being 
enhanced within the steps. Furthermore, it appears that 
different regions of the sample are at resonance at differ- 
ent values of the applied field, and that the sweep rate of 
the internal magnetic induction is spatially non-uniform. 
We demonstrate that the degree of non-uniformity can 
be manipulated by sweeping the magnetic field back and 
forth through part of the resonance. 

The local induction of single Mni2 crystals of typical 
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size 400 x 100 x 40 fim 3 was measured using an array 
of eleven Hall sensors 18 . The active layer of these sen- 
sors is a two dimensional electron gas formed at the in- 
terface of GaAs/AlGaAs heterostructure. The samples 
were mounted onto the surface of an array of 10 x 10 
fim sensors that were 10 /ira apart, with the easy-axis 
parallel to the z-direction and to the applied field H a , 
as shown in Fig. 1. The eleven sensors measure B x , 
the rc-component of the magnetic induction due to the 
magnetization of the crystal. Since the two dimensional 
electron gas resides only 0.1 \im below the surface, the 
induction measured by the sensors is practically equal to 
the induction at the crystal surface. As shown in Fig. 1, 
the center of the crystal was aligned with the array so 
that ten sensors probed the field along half the crystal 
length and the last sensor measured the field at the edge 
of the sample. 

Fig. 2a and 2b show two local hysteresis loops of 
B x (H a ), at T — 0.3K, measured simultaneously at the 
edge of the sample and close to its center. The curves 
were obtained starting from a fully magnetized state by 
sweeping H a from — 6T up to 6T and back. The loop 
measured locally at the edge of the sample (Fig. 2a) re- 
sembles previously measured magnetization loops. When 
increasing H a from zero to 6T, B x increases monotoni- 
cally toward its positive saturation value while display- 
ing a series of steep steps at roughly equal intervals of 
magnetic field, due to resonant tunneling of the magne- 
tization. 

In contrast, the magnetic induction measured close to 
the center of the sample (shown in Fig. 2b) displays 
non-monotonic behavior as a function of H a . The reso- 
nances are not manifested by a series of steps, but rather 
by a series of peaks and dips whose width is larger than 
that of the steps in Fig. 2a. Similar behavior is ob- 
served at other sensors located between the edge and the 
center of the sample. In general, B x is largest close to 
the edge and decreases towards the center of the sample. 
The non-monotonic dependence on H a , however, is most 
pronounced close to the center and gradually diminishes 
when proceeding towards the edge of the sample. The 
remarkable difference in the dependence of B x on H a ob- 
served by the different sensors contradict the assumption 
that the magnetization is uniform inside the sample. In 
general^, B x (f,H a ) = J dr'M(r',H a )f(f— r'), where 
M(r',H a ) is the magnetization at point r', that is as- 
sumed to be parallel to the easy axis, and 



f(f-P) = - 



3(z - z')(x - x') 
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If M = M (H a ) is independent of r', we get B x (f, H a ) = 
Mo(H a ) J dr'f(f— r'), so that the B x {H a ) curves mea- 
sured by all sensors should collapse onto a single curve 
when multiplied by a scaling factor. As the inner curves 
are non-monotonic while the edge curves are monotonic, 
this scaling cannot be realized, and the magnetization is 
therefore non-uniform. 
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FIG. 2: Local hysteresis loops of B x as a function of H a (T = 
0.3K) , measured simultaneously at two different positions on 
the sample: (a) close to the edge and, (b) close to the center 
of the sample. 



The non-uniformity of the magnetization is the source 
of the non-monotonicity displayed in Fig. 2b. Spins 
pointing in the same direction that are located on op- 
posite sides of a sensor make opposite contributions to 
B x . Due to these opposite contributions, the magnetic 
induction B x measured by a sensor located close to the 
center of the sample is very small when all the spins in 
its surrounding point in the same direction, but is en- 
hanced when the magnetization on the two sides is dif- 
ferent. This enhancement may bring B x to a value larger 
than its value when the sample becomes fully magnetized, 
as shown in Fig. 2b. Now, if the relaxation rates on the 
two sides of the sensor are different, the non-uniformity 
in the magnetization varies, and so does B x . The direc- 
tion in which B x varies with H a indicates the relative 
magnitudes of the relaxation rates on the two sides of 
the sensor. The situation is different for a sensor located 
near the edge of the sample, where all spins are on one 
side, and therefore their contributions to B x add. For 
this sensor, inhomogeneities in the magnetization only 
broaden the steps, but do not result in non-monotonic 
behavior. This qualitative picture is further discussed 
below in terms of an equation for the time evolution of 
the local magnetization. 

We now proceed to extract the spatial profile of the 
magnetization from the profile B x (z) measured by the 
sensors. Assuming that the sample has a perfect rectan- 
gular shape, and the magnetization is uniform along the 
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x and y directions, the induction B x can be written as 
B x (z) = - I ' dz'F{z-z')^p- (1) 

where 

cv t\ i w i + \A"i + rf i + ( z ~ z ') 2 
F(z - z ) = log v = - 

wi + \Jw\ + d 2 + (z~ z') 2 
w 2 + y/w% + d 2 + (z- z') 2 
w 2 + yjw\ +d 2 2 + {z- z') 2 

For the samples we used, d\ = 0.1 fim is the distance 
in the x direction between the sample surface and the 
2DEG, d 2 = 40 fim is the distance between the 2DEG 
and the sample's other surface, and w\ — —5 fim, w 2 = 
85 fim are the distances between the center of the sensors 
and the edges of the sample in the ^/-direction. 

F(z-z') is a short-range function, whose characteristic 
decay length is typically 10 [im for our geometry. Eq. (JTJ 
implies that the magnetic field B x measured by a sensor 
on the sample surface is proportional to the derivative 
in the vicinity of the sensor, and therefore Figs. 2a 
and 2b approximate the dependence of on H a at the 
edge of the sample and close to the center of the sample, 
respectively. 

Assuming that the magnetization varies smoothly over 
the scale of our sensors, we interpolate the measured 
B x (z) and numerically invert the kernel F(z — z') to ex- 
tract the magnetization profile M[z, H a ). Figure 3 shows 
the resulting magnetization as a function of position for 
various values of H a . Since our sensors cover only half the 
sample, we rely on the (nearly) symmetric shape of the 
samples and assume the magnetization to be symmetric 
with respect to the sample center. The non-uniformity of 
the magnetization and its evolution as a function of H a 
are clearly demonstrated. Focusing first on the central 
part of the sample, we note that at — 5.7T, where the sam- 
ple is presumably fully magnetized, the magnetization 
is approximately uniform. When increasing the field to 
3.6T the magnetization starts to develop a small amount 
of non-uniformity, manifested by the small variations of 
the magnetization from the average value. These varia- 
tions continue to develop and become more pronounced 
at 3.9T, when the magnetization in most of the sample 
shows a high relaxation rate. The non-uniformities then 
decay somewhat at 4.14T, where the sample is out of 
resonance, and become prominent once more at the next 
resonance, at 4.3T. When increasing the field further to 
5.7T the non- uniformities disappear and the magnetiza- 
tion becomes uniform again. Non-uniformities similar to 
those presented in Fig. 3 were observed in several sam- 
ples, where they are enhanced when most of the sample is 
at resonance, and become minor away from resonances. 

At the very edge of the sample the magnetization is 
expected to change abruptly from zero outside the sam- 
ple to a non-zero value inside the sample. We believe 
that the gradual variation observed at all field values re- 
sults from two major sources. First, the field B. x at the 
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FIG. 3: The calculated profile of the magnetization plotted as 
a function of the distance from the center of the sample, at six 
different values of the applied field. Note the enhanced non- 
uniformity at fields close to resonances, at 3.9T and 4.33T. 

edge of the sample varies significantly over the length of 
one sensor, and therefore the average field measured by 
the edge sensor may deviate appreciably from its local 
value. Second, the calculation of the magnetization is 
carried out assuming the sample has a perfect rectangu- 
lar shape. Deviations from this shape, for example edges 
that are not perpendicular to the plane of the sensors, 
could also contribute to the rounding of the data at the 
edges. Note, however, that these effects cannot introduce 
the observed non-uniformities in the center of the sample. 
We also note that the overall scale of the magnetization 
obtained from our calculation is not very different from 
that expected: at full magnetization, 4-7rAf is 0.12T for 
Mni2 crystals, while the maximum value shown in Fig. 
3 is 0.18T. 

The non-uniformity of the magnetization implies that 
B z , the internal magnetic field in the z direction, is spa- 
tially non-uniform. Consequently, different parts of the 
sample enter the resonance at different values of H a . We 
will now demonstrate that the internal sweep rate ^f- is 
spatially non-uniform as well, so that different parts of 
the sample spend different times at resonance. Under the 
assumption used to extract M(z), B z inside the sample 
can be written as 

B z (z) = H a + AnM(z) + f dr' dM{ r') z ' z ' (3) 

J oz' \r — r\ A 

Since all the terms and their time dependence are avail- 
able from our analysis, we can directly derive the local 
sweep rate ^gf-- 

This is shown in Fig. 4 as a function of position in 
the central part of the sample for different values of H a . 
Two major observations are of particular interest. The 
first is that at 3.9T and at 4.3T, when the sample is at 
resonance, the internal sweep rate is much higher than 



4 



in 



T3 




100 



FIG. 4: Calculated sweep rate of the internal field, ^§f-, 
plotted as a function of the distance from the center of the 
sample, for the same field values as in Fig. 3. Note that 
when most of the sample is at resonance, at 3.9T and 4.3T, 
the internal sweep rate ^f- is much higher than the external 
sweep rate, ^ £L . Moreover, at these fields is signifi- 

cantly non-uniform inside the sample. When the sample is 
out of resonance, at 3.6T and 4. IT, the internal sweep rate is 
approximately uniform inside the sample and is very close to 
the sweep rate of the external field. 



the nominal sweep rate of the external field. For exam- 
ple, is about 2.5 times larger than at 3.9T. In 
contrast, at 3.6T and 4. IT, when the sample is out of 
resonance, the sweep rates of B z and H a are approxi- 
mately equal. The second major observation is that at 
resonance fields, ^f- is highly non-uniform and can vary 
by as much as 20% over a distance of the order of 50 /im. 
This means not only that different regions enter the res- 
onance at different H a , but also that the times spent by 
local regions at resonance can be significantly different. 

As we show below, the observed spatial non- 
uniformities and their evolution with H a can be explained 
by a simple model describing the relaxation of the mag- 
netization along the hysteresis loop. Initially, at 6T, the 
magnetization is uniform inside the sample. However the 
relaxation rate is non-uniform due to the spatial variation 
of the internal magnetic field and due to molecular and 
structural disorder. Therefore, as we sweep the magnetic 
field to — 6T (in the opposite direction) non-uniformity 
in the magnetization begins to develop. Under certain 
conditions, and particularly within resonances, the non- 
uniformity in the magnetization grows with time, despite 
the decay in the magnetization itself. The magnetization 
is non-uniform even when the sample is out of resonance 
and the relaxation rates are negligible. We attribute the 
non-uniformity in the magnetization off resonance to the 



spatial variation of the sweep rate (Fig. 4) that causes 
different regions of the sample to spend different times 
at resonance. 

In our model the relaxation of the space- and time- 
dependent magnetization M(z,t) is governed by the 
equation 

_ = _ (M _M )r (4) 

where Mq is the equilibrium magnetization and T is the 
relaxation rate, which near a resonance has a general 
form 

[B z - B r y + a z 

where B r is the magnetic induction at the center of the 
resonance, and 2a is the resonance width. Molecular 
disorder introduces variations of r , B r and a between 
different molecules. However, since the distribution of 
these parameters does not depend on position, we use 
their average values. The non- uniformity of the magne- 
tization is characterized by its spatial derivative 0^- . We 
may use Eq. to analyze the time evolution of . To 
that end, we take the derivative of (4) with respect to z, 
and approximate B(z) rj H a + 47rM(z)(i.e., we neglect 
the third term in Eq. |3"|. which is small in the central 
part of the sample) . We obtain 



d dM 



= -[r + 47r(M 



, , . dT . dM 
M ^8B^ 



(6) 



dt dz 

Although the dependence of V on M(z) makes this 
equation a non-linear differential equation, the condi- 
tions under which non-uniformities in the magnetization 
grow with time are rather transparent. The question of 
whether the magnitude of Q£f increases or decreases with 
time is determined by the balance between the two terms 
on the right hand side of this equation. The first term 
acts to suppress the magnitude of i.e., the magni- 
tude of the non-uniformity. The effect of the second term, 
on the other hand, depends on the sign of (M — A1q)§^- 
When this sign is negative, this term tends to enhance 
the magnitude of the non-uniformity. If strong enough, 
it may overcome the effect of the first term. When that 
happens, the magnetization M(z) relaxes towards Mq, 
but its non-uniformity \ increases with time. 

From Eq. ©, it is easy to verify that as long as \M — 
M 1 



< — 



T + 4tt(M 
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(7) 
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and hence the non-uniformity is suppressed 
However, when \M — Mq\ > there is a 
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at which r + 47r(M-M )|g- < and the non-uniformity 
is enhanced with time. We consider, for concreteness, the 
case where H a is swept down from 6T to — 6T through 
one of the resonances. In this case Mo is negative, and 
M - M is positive. If M - M » ^r, the range © oc- 
cupies the first half of the resonance, where B z > B r . As 
H a is swept down, the non-uniformity grows during the 
first half of the resonance while B z > B r , and then is sup- 
pressed again when B z drops below B r . Such behavior is 
indeed observed in Fig. 3. In general, non- uniformities 
are always enhanced on approaching resonance, reach- 
ing maximum close to resonance and then suppressed af- 
ter the resonance. Our measurements show resonances 
whose half-widths are about 0.045T. However, the mea- 
sured half-widths are larger than the intrinsic half-width 
a, due to the effects of molecular disorder and hyperfine 
fields. The non-uniformities are found experimentally to 
evolve at the resonances where 4tt(M — Mo) > 0.06T, 
consistent with the model above. 

Eq. © and the discussion following it suggest that 
the non-uniformity in the magnetization may be en- 
hanced if the applied field H a is tuned to a value where 
T + 4n(M — M Q )^g is negative at least in some parts of 
the sample, and is kept there for some timeSI. Our mea- 
surements show that this is indeed the case. To identify 
the proper value of H a , we note again that B x is pro- 
portional to ^M- and therefore the observed slope j^f- 
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measures the local slope JjT^f • Thus, when H a is swept 
back and forth over a region where both and B x 
have the proper signs, this sweep results in an enhance- 
ment of the non-uniformity. Figure 5a shows the effect 
of such a back-and-forth sweep on the magnetic induc- 
tion B x measured by a sensor located close to the center 
of the sample. Two curves are presented: a reference 
curve in which H a was swept between — 6T and 6T at a 
constant sweep rate; and the 'back-and-forth sweep' (BF 
sweep) curve in which H a was swept up from — 6T to 3.9T 
and then was swept back and forth between 3.9T and 
3.69T. The reference curve shows a dip at 3.9T, signify- 
ing resonant relaxation of the magnetization. The range 
of the back-and-forth sweep was chosen to cover only half 
of this dip, at which both and B x were negative. 
During the back-and-forth sweep the magnetic induction 
\B X \ at H a = 3.9T increased to a value 3.5 times larger 
than its corresponding value at the reference curve. The 
strong increase in \B X \ above the saturation value signi- 
fies a strong enhancement of the non-uniformity inside 
the sample, since B x is proportional to the local spatial 
derivative of the magnetization, in the vicinity of 
the sensor. 

Figure 5b demonstrates the enhancement of the non- 
uniformity in the magnetization profile M(z) by the 
BF sweep. The lowest curve shows the magnetization 
at 3.68T just before sweeping the field back and forth 
over part of the resonance. At this point some non- 
uniformity has already developed. The three other curves 
correspond to three different measurements in which BF 
sweeping was carried out with three different amplitudes. 
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FIG. 5: Local hysteresis loops B x (H a ), measured close to the 
center of the sample. During the back and forth sweep the 
magnetic field B x at H a = 3.9T decreases to — 70G, which 
is 3.5 times larger than its corresponding value at the refer- 
ence curve (— 20G). The inset shows the evolution of B x as 
a function of time during the BF sweep, (b) Magnetization 
profiles M(z) showing the enhancement of the non- uniformity 
by the BF sweep. The lowest curve shows the magnetization 
at 3.68T before the BF sweep. The upper curves correspond 
to three different measurements in which the BF sweep was 
carried out with different amplitude of the BF cycle (0.08T, 
0.14T, 0.22T) 



In all measurements the cycle starts at the same field of 
3.7T, and each measurement spans a different part of 
the transition. The most significant enhancement of the 
non-uniformity is displayed for the smallest amplitude of 
0.08T, where the number of molecules that relaxed dur- 
ing the BF sweep strongly differs between the center and 
the edge of the sample. The effect is smaller for an ampli- 
tude of 0.14T, and yet smaller for an amplitude of 0.25T. 
While suggesting the possibility for magnetization non- 
uniformities that grow with time, Eqs. (JSJ) and © also 
indicate that eventually, as relaxation progresses and the 
magnitude of M - M gets smaller, T + 4tt(M - M )§^ 
must at some point become positive. Since Q£f is ap- 
proximately proportional to the magnetic induction B x 
measured by our sensors, this implies that when H a is 
swept back and forth for long enough, the magnitude of 
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B x must eventually decrease. This expectation is indeed 
borne out in the measurement. The inset of Fig. 5a 
presents the evolution of B x during the back-and-forth 
sweeping as a function of time. The measured B x ini- 
tially increases in magnitude, but after about 70 min- 
utes of sweeping, its magnitude starts to decrease. This 
change happens when the difference between the relax- 
ation rates on both sides of the sensor changes sign. The 
source of this change in relaxation rates may be the de- 
crease of the magnitude of M — Mo, or a change in the 
internal magnetic induction B z that changes both T and 
^g-. We are presently unable to determine the relative 
weight of these two sources. 

The saturation of the non- uniformity when M — Mq be- 
comes too small is plausibly also the source of the effect 
of the BF sweep on the following resonance. As shown 
in Fig. 5a, the BF sweeping of the field in the range 
between 3.69T and 3.9T changes the dip seen in the ref- 
erence curve at 4.3T into a step. Within the framework 
of Eq. ©, the non-monotonic behavior of B x around 
4.3T, shown in the reference curve, signifies an enhance- 
ment of the non-uniformity followed by its suppression. 
In contrast, the almost monotonic step observed in the 
BF curve indicates that as a consequence of the back- 



and-forth sweep, the resonance at 4.3T is traversed with 
r + A-rt{M — M )§jj remaining positive, such that no en- 
hancement of the non-uniformity takes place. 

To summarize, we used a set of Hall sensors to mea- 
sure the local magnetic response of crystals of the molec- 
ular magnet Mn 12 -acetate. Our measurements allowed 
us to determine the local magnetization in different re- 
gions of the crystal and its evolution as the crystal is 
driven through a hysteresis loop. We find significant non- 
uniformities, which are larger when resonant tunnelling 
takes place, but do not disappear when out of resonance. 
We explain how these non-uniformities result from the 
dipolar interaction between molecules. We show that the 
non-uniformity of the magnetization may be enhanced 
by sweeping the externally applied field back and forth 
through a properly chosen range. This method carries 
the potential for local manipulation of the magnetization 
profile inside the sample. 
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